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^ ■ Abstract 

■ In this work, we determine the short-distance coefficients for T inclusive decay into a charm 

\ pair through relative order within the framework of NRQCD factorization formula. The short- 

distance coefficient of the order-t;^ color-singlet NRQCD matrix element is obtained through match- 
ing the decay rate of bb(^S^^^) — )■ ccgg in full QCD to that in NRQCD. The double and single IR 
divergences appearing in the decay rate are exactly canceled through the next-to-next-to-leading- 
order renormalization of the operator 0{ S\ ) and the next-to-leading-order renormalization of the 
operators Oi^P^f^). To investigate the convergence of the relativistic expansion arising from the 
color-singlet contributions, we study the ratios of the order-f ^ and -v^ color-singlet short-distance 
coefficients to the leading-order one. Our results indicate that though the order-v^ color-singlet 
short-distance coefficient is quite large, the relativistic expansion for the color-singlet contributions 
in the process T —> cc -|- X is well convergent due to a small value of v. In addition, we extrapo- 
late the value of the mass ratio of the charm quark to the bottom quark, and find the relativistic 
corrections rise quickly with increase of the mass ratio. 

PACS numbers: 12.38. -t, 12.38. Bx, 12.39. St, 13.20. Gd 
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I. INTRODUCTION 

NRQCD factorization formula [l| provides a systematical approach to express the quarko- 
nium inclusive decay rate (cross section) as the sum of product of short- distance coefficients 
and the NRQCD matrix elements. The short-distance coefficients can be expanded as per- 
turbation series in coupling constant as at the scale of the heavy quark mass m. The 
long-distance matrix elements can be expressed in a definite way with the typical relative 
velocity v of the heavy quark in the quarkonium. 

The relativistic corrections to the quarkonium decay and production have been widely 
studied. In some processes, the next-to-leading-order (NLO) relativistic corrections are siz- 
able, some of which even surpass the leading-order (LO) contributions. For the quarkonium 
production, typical examples are the double charmonia J N) + rjc production at the B facto- 
ries 2|, |3| and the double rjc production through rji, decay ^. For the quarkonium annihilate 
decay, it happens in the J/ip inclusive decay and J/il){T) inclusive decay into a lepton pair 



or charm pair [5H7[. It may arouse worries about the relativistic expansion in the NRQCD 
approach. In Refs. j^, ^ , the authors resumed a collect of relativistic corrections to the pro- 
cess e^e~ J /ip + ?7^and found the relativistic expansion is of well convergence. Similarly, 
the authors of Ref. [4] studied the process 7]b — ric{mS)rjc{nS) and found the relativistic 
expansion is convergent very well, though the order-t;^ corrections are both negative and 
large. In Ref. 10|, the authors considered the order-w^ corrections to J/ip inclusive decay 
and found the contribution from the color-singlet matrix element in this order is not as 
large as the LO and NLO contributions. The order-^;^ corrections to this process were even 



calculated in Ref. where the convergence is furthermore confirmed. 

It has been shown that the order- 1;^ corrections to the process T — cc + X are extremely 
huge Ig]. The ratio of the short- distance coefficient of the order-f^ matrix element to that 
of the LO one approaches to —12. It seriously spoils the relativistic expansion. So it urges 
us to calculate the order-w^ corrections and investigate the convergence of the relativistic 
expansion for this process. Moreover, since the momenta of some gluons can be simulta- 
neously soft, there exist complicated IR divergences in the calculations. It is technically 
challenging to cancel the IR divergences and obtain the IR-independent short-distance co- 
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efficients through the color-octet mechanism ^, and it also provides an example to examine 
the NRQCD factorization formula. 

The remainder of this paper is organized as follows. In Sec. HIl we describe the NRQCD 
factorization formula for the T inclusive decay into a charm pair. In Sec. IIIH we list our 
definitions and present the techniques used to compute the decay rates. We elaborate the 
calculations on determining the short-distance coefficients at relative order v'^ in Sec. IIVI 
The QCD corrections to the color-octet operators are also calculated in this section. Sec. |V] 
is devoted to discussions and summaries. 



II. NRQCD FACTORIZATION FORMULA FOR T ^cc + X 

According to the NRQCD factorization formula, through relative order v^, the differential 
decay rate for T inclusive decay into a charm pair can be expressed as 



mP mP 



T 



, (o(^5P))t+ e '-^^{oi:^pf)h. (1) 

J=0,l,2 



where m indicates the mass of the bottom quark b, indicates the NRQCD 

-[c] 

'J 



operator with the spectroscopic state ^s+i^l^l ^[^\^ spin S, orbital angular momen- 
tum L, total angular momentum J, and color quantum number c, and {0{'^^~^^L^f))H = 
{H\0(^^~^^L^f)\H) indicates the NRQCD matrix element averaged over the spin states of 
H. The color quantum number c = 1 and 8 in the NRQCD operators denote the color-singlet 
and the color-octet respectively. The NRQCD operators in the factorization formula ([T]) are 



^ When we were doing this work, a work about ordcr-t;"* corrections to gluon fragmentation appeared at 
arXiv recently [l^ . In that work, the authors apphed the similar techniques to cancel IR divergences. 
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defined by 



22 



1 



^tya^ . xtr> 

1 ... . ^ 



^^'Tx ■ xU - t^) W + - 1^) Vx ■ X^(T^P 



(2a) 
(2b) 

(2c) 
(2d) 
(2e) 
(2f) 

(2g) 
(2h) 

(2i) 



where ib and y are Pauli spinor fields that annihilates the bottom quark and creates the 
bottom ant^uark. .espect.vel, denotes the PauM ,„at.x. S .ndicates the gauge-cova..,a,. 

derivative, and d = A — 2e represents the space- time dimensions. In (El), A^''B^^ and A^''B^^ 
indicate the antisymmetric tensor and the symmetric traceless tensor respectively, which are 
defined via [l^ 



2^ 



-{A'B' + A'B')-j 



(3a) 
(3b) 



In ([T]), we omit the term associated with the matrix element of the operator 



V^Vx -x^^ -gE + gE- ^)(T^ -gE + gE- 1S)ctx ■ X^(t^ 



(4) 



which is shown to be dependent and be a linear combination of the matrix elements of the 
operators fl2hl) and ([21]) accurate up to relative order 101 ] . 

In order to get the decay rate, both the NRQCD matrix elements and the short-distance 
coefficients in ([1]) should be determined. The NRQCD matrix elements have been studied 



through many nonperturbative approaches, such as lattice QCD [15[, nonrelativistic quark 
model 13], and fitting the experimental data [itI. \\\. 



The dimensional reg ularization in quarkonium calculations including definitions of operators was first 



given in Refs. [l3|, [14 \ 
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On the other hand, based on the factorization, the short-distance coefficients can be 
perturbatively determined through matching the decay rates of the relevant processes at 
parton level in full QCD to these in NRQCD. At the leading order in the coefficients 
dF{^S^^), (iF(^5'f'), and dF{'^PP) can be determined through the processes bb{^SQ^) — )■ ccg, 
bb{^Sf^) — cc, and bb{'^Pf^) — ccg respectively. The color-singlet short-distance coefficients 
can be determined through the process 66(^S'f^) — )■ ccgg. Among these short-distance co- 



efficients, -F(^S'|^') has been calculated up to the next-to-leadingorder in in Ref. 19 |. 



and F2{^sP) are obtained in Refs. [20|-l2^ and Refs. |6|,l7|, respectively. 
As mentioned in Ref. 10|, since there exists the relation {Q2(^SP))r = {Qi{^S^^^))r{l + 
to order v^, we can only determine the sum of the short-distance coefficients: ^ 

dFi'S?) = dF,{'S?) + dF.CS?). (5) 

As we will see, in order to get the IR-independent coefficient dF{^S^^), we are required to 
calculate dF{^S^i^) and dF(^Pf^), and make use of the color-octet mechanism to cancel the 
IR divergences appearing in the decay rate of the process bb{'^sf^) — ccgg. 

III. KINEMATIC DEFINITIONS AND PHASE-SPACE DECOMPOSITIONS 
A. Kinematic definitions 

We assign to the mass of charm quark, and assign pi and p2 to the momenta of 
the incoming bottom quark b and antiquark b. The momenta satisfy the on-shell relations: 
Pi = P2 = fTT-'^- Pi and P2 can be expressed as linear combinations of their total momentum 
P and half their relative momentum p: 

p,=P/2+p, p2 = P/2-p. (6) 

We assign Zi, I2 to the momenta of the final charm pair. Therefore the momentum of the 
virtual gluon yields to Q = + ^2- In addition, we take ki to the momentum of the gluon 
in the process bb — )■ ccg, and take ki, k2 to the momenta of the two gluons in the process 
bb — ccgg. 



In Refs. 



24 



-26[, the authors may provide a potential approach to distinguish dF^i^S^i^) and dF4{^S^i^) 
in ([5]), nevertheless, it is enough for us to determine dF(^S^^) in this work. 
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To facilitate the evaluation, we introduce a set of dimensionless variables 

_ 2h-P _ 2k2-P _4m| 

where M denotes the mass of the bottomonium. At the leading order in v, there is M = 
2m. All the Lorentz invariant kinematic quantities can be expressed in terms of these new 
variables. 



B. Phase-space decompositions 

In this subsection, we present the techniques for phase-space calculations of the relevant 
processes. We decompose each phase-space integral into two parts, which is proved to 
significantly simplify the calculations in the following section. 



1. 66(3 sfl) ^ ccgg 

The process 66(^5'^) — ?■ ccgg involves four-body phase-space integral, which can be ex- 
pressed as 

J y (27r)'^"i2A;? (27r)'^"i2A;0 (27r)'^-i2Z? (27r)'^-i2/o^ ^ ^ i 1 2; V ; 

As the treatment in Ref. j^, we decompose the space-space integral into 



dz 

^4 = /u</^3_l/U</^2-l, 



(9) 



where Jc?03_i and Jd(f)2-i are defined via 



On the other hand, we can also separate the squared amplitude of the process into two 
parts: the charm part and the bottom part, i.e., 

\M{^sf)\^ = L^'''^^''H\^\^s\:^), (11) 

pol.col 
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where polarizations and colors of the initial and final states are summed, a, b denote the 
color indices, and the charm part L^^^^^'^ is given by 

^iab),u ^ ^|^Tr[(K + m^h^ir, - m,)7l, (12) 
and the bottom part Hjj^i!'\^ S^i^) accounts for the remainder. According to the current 



conservation, we have 



ld<p,.,L'~^'^^'^ = ^ (- g^'^ + ^) X L, (13) 



where the Lorentz invariance L is explicitly calculated 



^ = ^(2 + -)Jl--. (14) 

As we will see, since L is a common factor in the whole calculations, we get (fT^ in 4- 
dimensions. As a consequence, the decay rate yields 

where a symmetry factor | is included to account for the indistinguishability of the two 
gluons. The second term in the parenthesis of (fT5|) does not contribute, due to the current 
conservation. 

It is accustomed to reduce Jd(j)3-i into the integration over two dimensionless variables: 



I' 



^3-1 = i287r3 r(2 - 2e) f^^'^'^A^^^ + X2 - 1 + ^)[(1 - Xi)(l - X2) - z] j. (16) 



where /, = (^e-^^)^ Here /i in represents the dimensional-regularization scale, and 
7E denotes the Euler constant. The boundaries of z,xi, and X2 are readily inferred 

1-xi- z<X2< ^ 0<xi<l-^, r<z<l. (17) 

1 — Xi 

It is convenient to make a further change of variables: 

This change of variables is particularly useful in extracting the relativistic corrections. 
Through this transformation, the boundaries of the new variables are simplified to 

0<a;<l, 0<y<l. (19) 



Finally, the decay rate reduces to 

r(^M ) = T^fe — Idzdxdy 



r(2-2e) 



1 - (1 - z)x 



(1 - z)-^'L X H('S\'^), (20) 



where 



2 (47r)4 1 - (1 - 2 



(21) 



Analogously, we separate the phase-space integral of this process into 

dz 

■'3-2 ^ I 2^ I "'/^2-l/ "(/^2-2, 



where Jc/02-i is given in ( llObl) and jd(j)2~2 indicates 

d'^-^Q d'^-^k 



(22) 



(23) 



' (27r)'^-i2go (27r)'^-i2fcO 
jd(f)2-2 can be straightforward integrated out as 

^ _ (1 - ^)^-^- r(i - e) 

^ ^ 8tt r(2 - 2e)r(l + e)'^' 

Similarly, we separate the squared amplitude into the charm part and the bottom part 
Huti^PP), and obtain the decay rate 



(24) 



r(i-6)/,e 

r(2 - 2e) 



Jdz{l-z)-^'L X HCpf^), 



where 



167r2 



(25) 



(26) 



3. bb{^sf^)^cc 



The phase-space integral of this process is quite simple. Therefore we directly present 
the decay rate 



r(35f ) = Jdz6{l - z)H('sf^) X L, 



(27) 



where H{^sf^) is defined via 



2 P2 



(28) 
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IV. DETERMINING THE SHORT-DISTANCE COEFFICIENTS 



In this section, we determine the differential short-distance coefficients of T include decay 

.[8 
J 



into a charm pair at relative order v\ dF{^Sf^), dF{^Pf), and dF{^S^l^) can be determined 



through calculating the decay rates of the perturbative processes hhi^Sf^) — )■ cc, hhi^Pf^) — )■ 
ccg, and 66('^S'f' ) — ccgg, respectively. 



A. S'-wave color-octet 

We first determine the differential short- distance coefficient dFi^sf^) through matching 
the decay rate of hbi^sf^) — )■ cc in full QCD to that in NRQCD. We also carry out the 
computations for the renormalization of the operator Oi^Sf^), which will produce mixing 
with the operator 0(^pf^) at the next-to-leading order in as and with the operator 0('^«S'[^') 
at the next-to-next-to- leading order in a^. Moreover we consider the renormalization of the 
operators Oi^Pf^)^ which will induce mixing with the operator Oi^S^^^). The renormalized 
operators will be utilized to cancel the IR divergences through the color-octet mechanism 
when we determine the differential short-distance coefficients dFi^Pf^) and dFi^S^^^). 



1. dF{^sf^) 

The corresponding factorization formula for hhi^sf^) — )■ cc is expressed as 

dV{'Sf) dF{^sf^) 

- ^^^^ 

Through calculations, we obtain the expression of H{^Sf^) defined in (^E^l as 

Hes?) = (30) 

Inserting fl5U]) into (^7^, we immediately deduce the differential decay rate in full QCD. 
By making use of (0(^5^ ))^^^3^[8]^ = 1, where the quark pair state is non-relativistically 
normalized, we readily get 

^£(;^ = ,(l_,)^xL. (31) 

dz 2 
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FIG. 1: The Feynman graphs for the NLO QCD corrections to the operator 0{^sf' ). We suppress 
the graphs which do not give rise to operator mixing. 



2. 



35p] ^ 3pj8] 



In this subsection, we consider the NLO QCD corrections to the operator 0{^Sf^). There 
are four diagrams illustrated in Fig. [1] ^. The vertex in the middle signifies the operator 
Oi^SP). Since the dia grams involve UV divergences, the operator needs to be renormalized. 
In this work, we uniformly use MS scheme to carry out renormalization. We express the 
renormalized operator as 



(32) 



where we truncate the expansion up to order a^, which is enough to current work. 6iO and 
62O are the corresponding NLO and the next-to-next-to-leading-order (NNLO) counterterms 
respectively. For convenience, we can expand an operator O in terms of a^: 



O = 0(0) + O(^) + + . . . . 

where the superscript '(n)' represents order-a" contribution. 

In the following calculations, we will also utilize the color decompositions [2 



(33) 



rparj-ib ^ rparpb 



4iV2 



N, 



■T" ® T" 



1 TV^ _ 2 

-1 ® 1 + ^, T" ® T". 



2Nr 



(34a) 



^ Here we consider only the diagrams which are relevant to our current work. Other diagrams will take effect 
when one considers the NLO QCD corrections to the short-distance coefficients. The similar calculations 
can also be found in Refs. 27 
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With the NRQCD Feynman rules 3l|], the color-octet contribution from Fig. [T]^a) reads 



4i(p ■ p' — p ■ kp' ■ k/k?') 
+ ie 



2m? Nc '^'^d-lj (27r)'^-i 

as T-®T-( I l\ , 

— p-p, (35) 



37rm2 Nc V^uv em, 

where = (47re~'^^)^. Since the coefficient is proportional to — the factor d coming 
from the loop integral in ( l35l) is replaced with 4 in the MS scheme. Moreover, we add the 
factor /e which is always associated with the MS scheme. 

It is not hard to find that Fig. [T](b-d) give the same contributions as Fig. [^a) up to a 
color factor. Summing all the contributions and employing ( 132|) . we get 



J 



3. 3p|'U3sl'l 



In the similar way, we get 



and 



OCPf'U = OCPf'r - ^ (i^I^OesI") + 0(al), (38) 
where Nj=l, '''^''j''''' , and for J = 0,1,2. In gUl and gg, we explicitly keep 

2 ■ - 

the factor d stemming from the S-wave extraction, i.e., p''p^ — t- ^^5*-'. The same treatment 
will be implemented in the following section ^. 



We thank G. T. Bodwin for communications on this point. 
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^. 3^[8]^35[1] 

We proceed to deal with the NNLO QCD corrections to the operator 0{^Sf^). The 
relevant Feynman graphs are illustrated in Fig. [2l There are totally eighteen Feynman 
graphs which do contribution. We are able to project out the color-singlet contribution by 
employing the color decompositions 

T'^T^T" (^T^T^T" = — l®l + others, (39a) 

27 ^ 

rpajbrpc ^ rj^bjarpc ^ ^ ^ ^ ^ + others, (39b) 

27 ^ ^ 

where we merely keep the color-singlet part. Through simple analysis, we find the color 
factors are ^ 1 ® 1 for the diagrams Fig. [2]^a,e-g,i,k,m,p,q) and — ^ 1 ® 1 for the diagrams 
Fig. [lb-d,h,j,l,n,o,r). 
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With the Feynman rules, the amphtude of Fig. [2]^a) reads 
{igsY f cfik f dH i 



'2a 



X 



po _ ko - /o _ + ^ep'O _ ^.0 _ /o _ (p'-fc-O^ + ,e 



4i(p p' -p kp' ■ k/k'^) Ai[{p - k) ■ {p' -k)-{p-k)- l{p' - k) ■ l/P] 

kP' + ie P + ie 

{d - 2)2 f d'^-^k f d'^"H [p ■ p'f 



4m4 (d-l)2 J {2txY~^ 7 (27r)'^-i |fc|3|Z|(|fc| + |/|)2' 
where the color factor is suppressed. 

Analogous to (HOil . we are able to get all others 

g^^ {d-2f f d'^-^k f d'^-H [p-p'f 



(40) 



he 



If 



Am^ {d - 1)2 J {271^-^ J (27r)^-i \k\Wi\k\ + |/|)' 
gt {d-2f f d'^-^k f d'^-H {p-p'f 



Am^ {d-lf 7 (27r)^-i 7 (27r)^-i |fc|3|Z|2(|fc| + |Z|) ^ 
gt (d - 2)2 f d'^-^k f d'^-H {p ■ p'f 



Am^d-iy J {2ttY-^ J (27r)^-i |fc|3|Z|3 ' 

h = Ik = In = hi 

h h Ig h Ip Iq h hi 

Ii = If- (41) 

Including the color factors, and summing over all these contributions, we get 

5 2g'^(d-2f f d'^-^k f d'^-H (p-p'^^ 



54m4(rf-l)2 ^ 7 (27r)'^-i 7 (27r)'^-i |fc|3|Z|3 



20«2 /e /. ' ' 



1 ® 1, (42) 



2437r2(c/- l)m4 \^euv 
where we use p'^ = p^. 

In addition, we also need to calculate the NLO QCD corrections to the counterterm 6iO. 
Employing fl37|) . we readily obtain 



gvrm^euv , 2437r2(d - 1) euv V^uv eiR/ 



J 

Combining ( H2l) and ( H3l) . we get 
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where the UV divergence can be canceled through the counterterm. Finally, we present the 
renormalized operator and the corresponding counterterm 

B. P-wave color-octet 

In this subsection, we determine the short- distance coefficients dF(^Pf^) through cal- 
culating the decay rates of the perturbative processes bb{^pf^) — )■ ccg. The factorization 
formulas for these processes are expressed as 

We can utilize ( l25l) to calculate the decay rates in full QCD. First, we need to obtain 
H{^pf^) defined in ( 1261) . To understand the IR structure and show the IR cancelation, here 
we separate Hi^pf) into two parts: H{^pf) = Hd(^PP) + Hsi^pf^). Hdi^pf) include the 
terms proportional to 1/(1 — z), which contribute the whole IR divergences to the decay rates 
from the region with the real gluon being soft. Hs{^pf^) take in charge of the remainder, 
which is absent of any singularity. Correspondingly, we use the subscripts d, s in both 
dr{^pf^) and dF{^pf^) to denote the contributions from the two parts. The advantage 
of this classification will be recognized when we determine the short-distance coefficient 
dFi^sP). 



1. H.epp) 

We directly present the expressions of Hd{^PP) 

H r3p[8li - 5ag(l - e)p^ 

> - 3m4(l-z)(3-26)' ^^^^ 

which are the same for J = 0, 1, 2. 

It is useful to make the following expansion 

1 (1-^)"^%/ X 1 ^ /ln(l-^) 



S{l-z)+ —-] -26^—^ +Oie% (48) 



14 



where the plus functions are defined via 



dz{a)+f{z) 



dza 



(49) 



Inserting (1471) into (125|) and employing (l48l) . we are able to get the differential decay rates 
dVdi^Pf^), which embrace IR divergences. Incorporating the factorization formulas fH6l) and 
the expressions (13T]) and fl36al) . we find the IR divergences appearing on the left hand side 
of (H6l) are exactly canceled by the renormalized S'-wave color-octet matrix element on the 
right hand side. It renders the short-distance coefficients dFi^Pf"^) free of any singularity: 



dz 



7r2 28 



^(l^ll\+^-2 ln(l - r) + (^2 In^fl - r) - — ln(l - r) 
.eiRV fej 3 V 3 



5(1 - z) 



1 - z 



lOe 



V 1 - ^ 



(50) 



In (ISO]) , we keep the terms linearly dependent on e, which will induce finite contributions to 
the short-distance coefficient dFi^S^^). 



2. Hs{^Pf) 



We present the expressions of Hsi^pf^) as 



Hs{'P^' 



HsCP?^) 



HsCp?] 



5a2(5-^)(l _e)p2 

12m4(3 - 2e) ' 
5al[2 + z - 2{1 + z)e]p'^ 
6m4(3 - 2e) ^ 
5a^[4 2 - (1 2;)e]p2 



(51a) 
(51b) 
(51c) 



6m4(15- 16e + 4e2) 

The corresponding differential decay rates can be achieved by multiplying a constant factor: 



dz 
dz 
dz 



5(5 - z)a'Me 
36 



l-2eln(l-^) + y 



" ■' 2 + 2 - 2e(^ 2) ln(l - 2) + ^ -e 



18 
18 



z + A- 2e(z + 4) ln(l - z) + —(31z + 169)e 

15 



(52a) 
(52b) 
(52c) 
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C. 5- wave color-singlet 



In this section, we determine the short-distance coefficient dF(^S^i^) through calculating 
the decay rate of the process bb(^S^i^) — t- ccgg. The corresponding factorization formula is 
expressed as 



mP rrr 



-^^^^^1^^(0(35^)), + $:^^^^ (53) 

J 

where the subscript 'if' in the matrix elements represents 66(^5'^). As mentioned in Sec. [Tll 
the two matrix elements in the second line of (|53l) are equal at relative order v^. Therefore, 
we will determine the combined short- distance coefficient dFi^S^^^ defined in ^ in this 
subsection. 

There are six diagrams for this process. The formula for the decay rate is given in (!20l) . 
Firstly, we need to calculate B.{^S^^\ which is defined in fl2T]) . To separate the relativistic 
corrections , we expand E{^S\^) in powers of p : 

H{^S\^) = if(o)(35W) + i/(2)(35W)^ + H^^\^SP)^ + 0(p6), (54) 

where the ffist two orders have been considered in Ref. j^. Our remainder task is to calcu- 
late H^^^'^sf^) and the corresponding decay rate. At order v^, the decay rate involves IR 
divergences. Analogous to the treatment in the previous subsection, we separate 
into three parts: H^^^^sP) = H^f\'sP) + Hi^\'sP) + H^^\''S^^^). //^'^(^ 5™) proportional 
to 1/(1 — 2;) contributes the whole IR divergences to the decay rate in the region where the 
two real gluons are simultaneous soft. In the following, we will demonstrate that the IR 
divergences can be thoroughly canceled by the renormalized S'-wave color-octet matrix ele- 
ment fl45ap with the short-distance coefficient flHT]) . together with the renormalized P-wave 
matrix elements fl38|) with the short-distance coefficients dFd{^pf^) given in 
contributes the whole IR divergence to the decay rate in the region where only one of the 



^ Since in the region z — s- r. the decay rate wih develop a logarithmic dependence on r, i.e., Inr, which is 
sensitive to the value of r, we will not expand the r appearing in L. For further explanations, we refer 
the reader to Ref. 
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real gluons is soft, as a result, it should be proportional to either 1/x or 1/(1 — x). We 
will show that the IR divergence can be thoroughly canceled by the renormalized P-wave 
matrix elements ( 13 8 p with short-distance coefficients dFsi^Pf^) given in f l52p . 
takes in charge of the remainder and therefore corresponds to a finite contribution to the 
decay rate. Similarly, we use the subscripts d, s, r in both rfr(34^]) and dF{^S^^^) to denote 
the contributions from the three parts. 



1. Hfi^S^?) 



In this subsection, we calculate the decay rate and short- distance coefficient related to 
H^f^i^S^i^). First, we present the expression 



3(22/2 - 2y + 1) + (8y2 - By + l)e - 8y(l - y)e' 



.(55) 



243m27rx(l -x)(l - z) 

This term is proportional to 1/(1 — 2;). In addition, we notice it also contains the factor 
l/x(l — a;). Therefore, we expect it will contribute a double IR pole to the decay rate at 
the endpoints of z and x. 
Employing the expansion 



:i-r) 



-4e 



l+4e 



4e 



-6{l-z) + 



4e 



ln(l 



1 



0{e% (56) 



and integrating out the variables x, we are able to obtain the differential decay rate 



dT,{^S^^\ 
dz 



20agL/,V 
2437rm6 

x5{l 



1 ^ 4 -41n(l 



z - 



- + 2\\iz 

e 



+ 81n2(l -r) - 161n(l 
ln(l - z) 
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+ 16 



35 

y 

(57) 



Inserting (157|) into the factorization formula (!53l) . we see the IR divergences on the 
left hand side are exactly canceled by the S-wave color-octet contribution, together with 
the renormalized P-wave color-octet matrix elements with the short-distance coefficients 



dFcti^Pj) on the right hand side. Straightforwardly, we get 



3pj8h 

dF,{^sl^) 
dz 



lOa^L 
7297r 
27r2 



In' 



2/^ 



10 - 121n(l 



5{l 



M2 



z - 



-^lnA_ + 41n^(l 



.1 1 20 

41n^ + 21n. + - 



ln(l - z) 



, 25 

r) H 

' 9 



, (5^ 



where indicates the factorization scale. 
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2. Hf\^sf) 



The expression oiEf^{^S\^) reads 



20af r-8?/2 + 8?/-^-3 2(z - 5)(?/2 - + ^ - 1 



1 



-e — 



1 — X 



X X (^1 — y + i/z)^ 

12(1 - zfy^ + 12(2 - 5)(2 - 1)V - 6(2^ - 11^^ + 31^ - 21)?/^ 
-24(^2 - 5z + 6)y^ + [z^ + + 9^^ - TJz + 96)2/^ + 2(2^ - 2z^ - 3z - 18)y 

+ 3;z + 6 + e(^16(l - 2)V + 16(^ - 5)(z - 1) V - 2(2 - 1) 
X (^2 - 342 + 81)?/^ - 8(^2 - 142 + 21)y3 - (2^ - 2^ - 9^2 - 32 - 92)i/2 



-2(32^ + 7^ + 12)1/ -2^ + 2 + 2 



(59) 



This term is proportional to either 1/x or 1/(1 — x\ therefore, it contributes single poles 
to the decay rate at the endpoints of x, i.e., x — > and x — )■ 1. After integrating out the 
variables x, y, we have 



dz 



7297rm6 [ e 3(1 - 2)^ 
-3(32^ - iz^ + 42^ + 122^ - 92 + 6) In z 



2l2^ - 43/ - lOlz^ + 291;^^ - 248;^ + 80 



- 4(3z + 5) ln(l - 2) 



(60) 



It is not hard to find that the IR divergence in (!60l) is exactly canceled by the renormalized 
P-wave color-octet matrix elements with the short- distance coefficients (iF,(^p'^^ 
expected. We present the final result 



^ J as we 



dz 



1 



1^1(3. + 5) ln^ + 

2 1 877r\^^^ + ^^^''M2 + 3(l-2)4 

+3(32^ - 4/ + Az^ + 122^ - 92 + 6) In 2 



62^ - 16z^ - 59z^ + 153z^ - 131;z + 47 



2(32 + 5) ln(l -2) 



(61) 



3. Hi^''\^S?) 

Finally, we deal with Hr'^\^S^i^), which will produce a finite contribution to the decay 
rate. Since the expression of Hr'^\^S^^^) is both tedious and cumbersome, here we merely 
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present the differential short- distance coefficient as 



dz 



2624407r(l - zY 
-56942 + 1609)tan~^ 

-72962=^ - 4762^ 



15tan ^ 



1 - z 



eCSe^^ - 279^'' - 25482'^ + 6718^^ 



1 - z 



+ ^z{l - 2)(13908 - 382152 + 33375^ 



6(2992^ + 226402^ - 278482^ + 162082=^ - 20152 + 700) 
X In^ - (1 - 2)(68916z^ - 256073^^ + 696819^^ - 669618^ + 229580)|. (62) 

D. Summarizing the differential short-distance coefficients 

We summarize the differential short-distance coefficients calculated in the above: 

(63a) 
, (63b) 

(63c) 



dF{^sf^) 
dz 


= ^"^^5(1 
2 ^ 






dF{^pf) 


-balL / 




5 


dz 


" 18 1 


+ 3- 


dF{^sf) 
dz 


_ dF,{^s[^ 
dz 


) dF, 
- + — 


dz 



21n(l -r 



5{l-z)-2 



1 - 2 



A, 



dz 

where Aj are — 2 + 2;, and for J = 0, 1, 2 respectively. 



V. DISCUSSIONS AND SUMMARIES 
A. Discussions 

Applying the formulas (!63|) for the differential short- distance coefficients obtained in the 
last section, we now make some discussions. 



1. dF(^sP)/dz in the limit of z ^ 



We first discuss the short-distance coefficient dF{^S^i^)/dz in the limit of 2 — 0. It is not 
hard to derive 

1 dFi'S?} 



L dz 



35a3 70^3 „2 1609^3 3913af 
Inz In- ' 



(64) 



21877r""' 2 1 877r M2 11664 43747r 
In (IMj) . we notice that the limitation bears the logarithmic divergence In 2. It is no surprise, 
owing to the fact we actually do not regularize the singularity when z approaches to 0. 
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FIG. 3: Ratios of the differential short-distance coefficients. The blue dashed line represents the 
distribution ti{z), while the red dot-dashed line represents the distribution t2{z). We specify the 
factorization scale fj, = M in t2{z). 

Moreover, we find that the coefficient of In z equals exactly to the corresponding coefficient 
of IR divergence in the decay rate of the process 66(^5'f' ) — )■ "ig up to a factor x |, as 
our expectation (The constant factor originates in from (fT5|) .). On the other hand, when 
either of the two real gluons becomes soft, there exists IR divergence which is regularized 
in dimensional regularization and canceled by the renormalized P-wave color-octet matrix 

2 

elements. The term proportional to In^pj is related to the divergence. We are delight with 

2 

that the coefficient of In^^ is exactly double that of In z. 
2. Color-singlet differential short- distance coefficients 

In (1^ . the two types of plus functions and C^^^Z^^ )+ diverge as 2 — t- 1. Since 

these singularities arise when the momenta of the real gluons in the final states go to 0, 
the distributions are actually unreliable in this region. Nevertheless, the singularities in the 
distributions are smeared when one integrates out z and so the integrated short- distance 
coefficients are well behaved. In order to investigate the dependence of the relativistic 
corrections on the virtuality of the intermediate gluon, it is intriguing and enlightening 
to study two ratios: ti{z) = dF2Csf^)/dFi{^sP) and t2{z) = dF{^sf^)/dFi{^sP), where 
dFii^S?) and dF^i^'S?) are 

see it clearly, we plot Fig. [3] to illustrate the two ratios. In the figure, we observe that ti{z) 
which reflects the NLO relativistic corrections is negative and its magnitude rises quickly 
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TABLE I: Ratios of the short-distance coefficients for T ccgg. The charm mass is selected to 



that of the D meson [g, 



221 ]. The masses of the D meson and bottomonia are taken from Ref. 



r 


F2 






T(15) 1.56 X 


10-1 


-12.4 


19.4 + 0.61n2(i^) + 0.91n(^) 


T(25) 1.39 X 


10-1 


-11.7 


18.5 + 0.5hi2(^) + 0.6 In(^) 


T(35) 1.30 X 


10-1 


-11.4 


18.0 + 0.51n2(^)+0.41n(^) 



with increase of the virtuality of the immediate gluon. However t2{z) which reflects the 
order-f ^ relativistic corrections is positive in small values of z and turns to negative in large 
values. The magnitude of t2{z) is sizable in most values of z. 



3. Integrated color-singlet short- distance coefficients 

Finally, we integrate out the variable z and investigate the integrated short- distance 
coefficients. In Table. [U we list the ratios of the order-w^ and the order-?;"^ color-singlet 
short-distance coefficients to the LO one for the processes of T{nS) inclusive decay into 
a charm pair. We learn from the table that the ratio F(^S^i^) / Fii^S^i^) is both positive 
and sizable, nevertheless, the relativistic expansion is of well convergence due to a small 
value of V (e.g., ~ 0.1.). The situation is quite similar to the case for the process 
J /ip ^ ggg [lOj]. To further study the relation between the relativistic corrections and r, we 
extrapolate the value of r and investigate the ratios: gi{r) = ^2(^5^' )(r)/Fi(^S'|^')(r) and 
g2{r) = F{^S^^){r)/Fi{^sf'^){r). We illustrate the two functions in Fig. HI From the figure, 
we find the NNLO relativistic corrections become more and more important with r increase. 



B. Summaries 



In this work, we determine the short- distance coefficients within the framework of 
NRQCD factorization formula for T inclusive decay into a charm pair through relative or- 
der v^. The order-w^ color-singlet differential short-distance coefficient dF{^sf^^) is obtained 
through matching the decay rate of bb(^S^i^) ccgg in full QCD to that in NRQCD. The 
double and single IR divergences appearing in the decay rate are exactly canceled through 
the NNLO renormalization of the operator 0{^Sf'^) and the NLO renormalization of the 
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r 

FIG. 4: Ratios of the short-distance coefficients as functions of r. The blue dashed hne represents 
the distribution gi{r), while the red dot-dashed line represents the distribution g2{r). We specify 
the factorization scale fj, = M in g2{r)- 

operators Oi^P^f^). To investigate the magnitude of the relativistic corrections and the 
convergence of the relativistic expansion, we show both the ratios of the differential short- 
distance coefficients ti{z) and the ratios of the short-distance coefficients gi{r). Our results 
indicate that though F(^S'P) is quite large, the relativistic expansion from the color-singlet 
contributions in the process T — ?■ cc + X are well convergent due to a small value of v. In 
addition, we extrapolate gi{r) to a large range of r, and find the relativistic corrections rise 
quickly with increase of r. 
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